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It is now well established that Zn doping of high-Tb cuprates reduces their Tc and triggers the 
appearence of a spin glass phase. In this context, we have solved exactly the problem of N non 
magnetic impurities in the staggered flux phase of the Heisenberg model which we assume to be a 
good mean-field approximation for the spin-gap phase of the cuprates. In this model, the quasi- 
particule spectrum has four nodes on the Fermi surface, and linearization of the spectrum in the 
neighbourhood of these nodes leads to a system of 2D Dirac fermions. In the presence of a macro- 
scopic number of (non magnetic) impurities, the problem has a characteristic logarithmic structure 
that renders ineffective the usual perturbative expansions. We have used this logarithmic structure 
to calculate an exact solution. For a concentration ni of impurities in the unitary scattering limit, 
the additional density of states is found to be proportional to ni/(w ln 2 (\w\/D)) (where D is the 
infrared cut-off of the linearized spectrum) in analogy with the ID case of doped spin-Peierls and 
two-leg ladders compounds. We argue that the system exhibits a quasi long-range order at T = 
with instantaneous spin-spin correlations decreasing as ni / In 4 (ni / Rij) for large distances Rij be- 
tween impurity sites. We predict enhanced low energy fluctuations and compare these results to 
NMR and inelastic neutron scattering experiments in the high-Tc cuprates. 
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It is now established that the normal phase of the un- 
derdoped cuprates possesses a magnetic gap with the 
same anisotropy as the d-wave superconducting gap. The 
essential physics of the high-Tc* superconductors can be 
captured by focusing on the Cu02 planes. A microscopic 
starting point of theoretical analysis is the so-called t-J 
model. In a slave-boson mean-field formulation, the sys- 
tem undergoes spin and charge separation: an electron 
in these highly correlated materials is a composite object 
made of a spin 1/2 neutral fermion (spinon) and a spinless 
charged boson (holon). The gap in the magnetic excita- 
tions of the normal state can be viewed as singlet for- 
mation between pairs of neutral fermions in the absence 
of coherence between the holons. We assume that this 
mean-field picture captures essentially the initial idea of 
Anderson jlj] of a Resonance Valence Bond ground state 
for the normal phase of the cuprates. 

Substitution of Cu ions in the conduction planes of 
high-Tc cuprates by different non magnetic ions presents 
an important experimental tool for the study of the 
metallic state. Unusual effects has been revealed espe- 
cially when the materials were doped with Zn. The va- 
lency of Zn is Zn 2+ (d 10 ) and compared with the Cu 2+ 
case, one electron is trapped by an additional positive 
charge of the nucleus, forming a singlet at the Zn site. 
In the spin-gap phase, it is experimentally found by 
NMR p pj experiments that a local magnetic moment 
of spin 1 [2 appears on Cu sites neighbouring the Zn im- 
purity. 

In high-Tc cuprates the persistance of AF fluctuations 
in the metallic state is probably one of the most striking 
feature. Inelastic Neutron Scattering have established 
previously the existence of an energy gap in the imagi- 
nary part of the dynamical susceptibility in the normal 
and superconducting phases of the pure compound ||] 
(without Zn). Substitution of Zn in Cu02 planes shows 



a striking transfer of spectral weight from high to low 
energies, partially filling the spin-gap This is the 

signature of strong enhancement of AF correlations in 
the spin-gap phase. 

Earlier work has treated the effect of a single non 
magnetic impurity in the 7r-flux state || and the pair- 
wise interaction between them [ju^]. It was found that 
each impurity creates a bound state in the pseudogap 
at lo = 0. Interaction between a pair of impurities 
leads to a level splitting between these states given by 
A(i?) ~ l/(R\n(R)). 

Here we solve the problem of N non magnetic impuri- 
ties in the staggered flux phase of the Heisenberg model. 

In the mean-field of the t — J model at half filling, the 
7r-fiux phase is a system of 2 dimensional Dirac Fermions 
with nodes located at the for points (±7r/2, ±7t/2). The 
impurities are treated as repulsive scalar potentials ran- 
domly distributed on the lattice. 

The key idea of this solution is to write the T-matrix 
equation directly for N impurities and then to factor- 
ize the leading divergence. As the impurities are treated 
as scalar scattering centers, equations of motion can be 
closed algebraically by a N x N matrix M which carac- 
terize the problem. 

Here My oc 4/ (irD 2 )iu> n Ko (Rij\ui n \/D) if i and j be- 
long to the same sublattice and 

My oc 4/ '(nD 2 )iw n Ki (Rij\ui n \/D), if i and j belong 
to different sublattices. Now T = — VqM^ 1 and in 
the unitary limit, additional density of states can be 

casted into the form Sp(uj) = —XjixlmTr 

-Vo/irlmTr (m~ 2 8M 2 /du\. 

All the configurations of impurities produce a loga- 
rithmic divergence in M 2 that we factorize : M 2 oc 
In \(J n Rij/D\ S, where S is a N x N matrix satisfying 
Sij = 1 if < D/\oj n \ and elsewhere. We argue that 
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the leading part of the density of states is due to the 
logarithmic divergence and the effect of S is negligeable. 
The value of S above is an approximation by excess of 
the rest. 

We take = (\u n \/D) Uij. For all j inside a circle 
of radius D/\u n \ around the point i, Uij is a random 
configuration of ±1 so that J^j Uij ~ D/\u n \. In addition 
all the points j situated in the external boundary of this 
circle have Uij = Elsewhere Uij — 0. The main 

difficulty in the inversion of S is that two circles of radius 
D/\u n \ centered around two points i and j very close to 
each other will overlap, leading to the same number of 
non zero coefficients in the lines i and j of S. In order 
to differentiate the sums ^ fc SikS^ 1 and ^2 k SikS^ 1 we 
have thus used the external boundary of the circle to 
compensate its volume in S* -1 . 

With this inverse we show that the effect of S is neg- 
ligeable in the density of states. The result is exact in 
the limit of low frequencies and in the unitary limit. We 
find 
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For one impurity alone P]lC|] , a S- like bound state is cre- 
ated at u = 0. If the impurities were totally uncorrelated, 
we would find N 5-functions at u> = 0. Here we see clearly 
the overlap of impurity states which leads to a broaden- 
ing of the 5 functions. It is interesting to compare this 
density of states with the one found for one dimensional 
spin-gap systems such as spin-Peierls and two-leg ladders 
systems where Sp(uj) ~ 1/ (\lj\ In (|w|/D)). 

In order to adress the question of long-range order in 
this system, we compute the instantaneous sorrelation 
function (S+S~). We find that at T = 0, 



(In 2 imRij) + (tt/2) 2 



(2) 



for impurities in different sublattices. N non magnetic 
impurities randomly distributed in a 7r-flux phase sta- 
bilize a quasi long range staggered order. Physically 
we understand what happens by introducing impurities 
one after another. For one impurity, a local moment 
with S = 1/2 is created in the vicinity Two im- 

purities interact via an effective Heisenberg exchange 
H = J(R)S X .S 2 , with J(R) ~ \/{R\n(R)) if they are 
located on different sublattices. However two impurities 
located on the same sublattice don't interact. This com- 
mensurability effect with sublattices gives the intuition 
that no frustration will be introduced in the model when 
a macroscopic number of impurities interact with each 
other. This effect is the origin of the quasi long range 
order that we find. 

The instantaneous correlations on different sublattices 
are proportional to the density of impurities because 
a pair of impurities interact via an effective exchange 
J(R) ~ I/(i?ln(i?)). We suggest that the logarithmic 



decay in (g) is due to quantum fluctuations at zero tem- 
perature even though the precise power in the logarithm 
may depend on our approximation. 

What would be observed if a Neutron Scattering ex- 
periment would be performed on this system? We 
have computed here the structure factor S(Q, 57) = 
iM^Sg x"(Q) integrated around (7T,7r). Using the 
same tricks as before for the instantaneous correlation 
function, we show that at low frequencies S(Q, f2) ~ 
Uij (f21n 4 |fi/D|). Even though we don't have true long 
range order in the system (which would have led to 
S(Q, 51) = 6 (CI) ), Neutron Scattering experiments would 
see a divergence of the integrated staggered structure fac- 
tor at low frequencies, signaling a strong enhancement of 
antiferromagnetic correlations. 

The situation will change for experimental systems 
with finite oxygen doping. If we take as a theoretical 
starting point the d-RVB ground state, the effect of oxy- 
gen doping will be to move the nodes away from the four 
points (±7r/2, ±7r/2) in the canonical spectrum. Exci- 
tations from node to node will then slightly differ from 
the commensurated vectors (n,ir) or (tt, 0), (0, n). We 
believe that instead of quasi-long range AF order, this 
mechanism may give rise to an incommensurate quasi- 
long range order. Furthermore, it is likely that the finite 
frequency response integrated over g-space will not be 
that different from that we find, which may explain the 
enhancement of AF correlations observed in INS. 
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